NMR experiment factors numbers with Gauss sums 
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According to the legendQ Karl Friedrich Gauss as a child was given by his teacher the problem of adding all 
integers up to a given number m. His elegant solution rearranges the terms of the series in a convenient way arriving 
at the result 
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which depends quadratically on m. In his adult life Gauss revisited square dependencies in the context of number 
theory [2J when he calculated the sum of quadratic phase factors which today carry the name Gauss sums [3. In 
this contribution we present the first implementation of a factorization algorithm [J, |5j based on Gauss sumsjg] by 
factorizing the number N = 157573 using NMR techniques. It is intriguing that the physical realization relies on the 
sum formula Eq. Q . 

Factorization of numbers into their prime factors is considered a hard non-polynomial problem^ for classical 
computers. It was Shor Q who proposed a quantum algorithm which can solve the problem on a quantum computer^ 
with a tremendous speedup as compared to a classical computer. Vandersypen et al.[lf1] demonstrated in an NMR 
implementation of the Shor algorithm with seven qubits the factorization of the number 15 = 3 • 5. 

Gauss sums are the discrete version of the Fresnel integrals^]] familiar from classical optics effects such as diffraction 
from a wedge or the Feynman formulation of quantum mechanics. Thus they are a sum rather than an integral over 
quadratic phase factors. Due to this discreteness Gauss sums have interesting periodicity properties which manifest 
themselves in the Talbot effect ^3]> fractional revivals or curlicues 0]. -Due to this periodicity property they not 
only play an important role in number theory but are also an ideal tool to factor numbers. Indeed several such 
Gauss sum factorization algorithms [J, |a, Ujj UJl UJ1 Oil have been put forward. All these schemes capitalize on the 
periodicity properties but differ in their method of implementation ranging from a AT-slit Young interferometer^^ 
via molecules 0| and wave packet dynamics ^3 to chirped laser pulses interacting with atoms 0, ]jj 0|. However, so 
far no experimental demonstration of this approach has been provided 

In the present paper we use NMR techniques to realize a rather special Gauss sum which brings out factors with 
a remarkable contrast even when only a few terms appear in the sum. We start by briefly summarizing the essential 
ingredients of our scheme. Since we do not employ entanglement yet, we still need exponential resources. Nevertheless, 
the quasi-random interference |2d| contained in the Gauss sums allows us to work with only a few terms in the sum. In 
the experiment this advantageous feature translates into the need for a few pulses within a decay time of the system 
only. We then turn to the experiment and demonstrate the power of the scheme for a number with six digits, but 
emphasize that extensions to much larger numbers are within reach. 

Our experiment implements the Gauss sum|21| 
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with M + 1 terms and N is the number to be factored. The argument £ with 1 < £ < \ZjV scans through all integers 
between 1 and for possible factors. 

The capability of the Gauss sum, Eq. (J2J, to factor numbers stems from the fact that for an integer factor q of N 
with N = q ■ r all phases in A^ are integer multiples of 2n. Consequently the terms add up constructively and yield 



A^f\q) = 1 as shown in Fig. by black dots. When £ is not a factor the quadratic phases oscillate rapidly with 



N 

and A^ 1 ^ takes on small values. 
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FIG. 1: Factorization interference pattern for TV = 157573 = 13 • 17 • 23 • 31 obtained from the Gauss sum A 
and from the Gauss sum corresponding to damping (circles). In both cases already M + 1 = 11 terms are sufficient to 

clearly discriminate the factors from non-factors. Note the different scales for A^ 1 ' 1 and CL . 



In this process of destructive interference the truncation parameter M plays a crucial role. Indeed, already a few 
terms in the sum are sufficient to discriminate factors from non-factors in the signal |^4^^(€)|. In order to analyze 
the dependence of this surprising feature on M we define the contrast V = (l~ a )/(l + a ) °f the factorization 
interference pattern in complete analogy to classical optics where a = J2e ,< Li (^')\/ n o denotes the average value 

of the sum at the non-factors upto y/~N. Indeed, no is the closest integer to y/~N and the summation runs over all 
arguments I' which are not factors of TV. In Fig. [2] we show V for different values of TV. Already a relatively small 
number of terms M results in good contrast of the signal. 
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FIG. 2: The contrast V of the factorization interference pattern of Fig. Q as a function of the number M of terms in the Gauss 

(M 
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sum A^' n for TV = 157573 (circles) and TV = 4683359 (black dots) 



Next we address the complexity of our factorization scheme. To gain information on the factors of TV we have 
to measure the signal l-A^^^)! for arguments I belonging to the interval [0, n a ]. Since at most \fN data-points 
{I, |^4^ W|} have to be acquired, we estimate the required resources as \/TV = exp[L/2] where L — log TV is the 
number of digits of TV. Although our scheme scales exponentially we profit from the small number of terms M 
necessary to distill the factors. 

We now turn to our experimental realization of a Gauss sum. In the original proposal^ the Gauss sum arises 
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through the time evolution of a two-level atom whose transition frequency increases linearly in time and which is 
driven by a train of laser pulses. In our experiment we subject an ensemble of spins I = 1/2 to a specific sequence of 
RF pulses as shown in Fig. |3J They lead to a sequence of signals which when summed represent a Gauss sum closely 
related to Eq. J2J. 
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FIG. 3: NMR implementation of our factorization scheme using Gauss sums. The j/-pulse which prepares the initial density 
matrix p ln is followed after a time r by a pulse which imprints the phase cf>o on the spins. This pulse is the first of a sequence 
of M + 1 pulses each imprinting a different phase 4> m . At times r m + r = (2m + l)r + r we measure the echo, that is the 
polarization s m in :r-direction and sum the echos s m over all m to obtain S^ M ' > . In order to implement a Gauss sum the phases 
4> m need to be proportional to the number N to be factored and have to increase linearly as a function of m since the spin 
dynamics expressed by s m depends on the sum over all phases of the previous pulses. 



In our experiment we use H2O as an ensemble of protons with nuclear spin 1/2 in Boltzmann equilibrium at room 
temperature described by a density operator pg = t/2 — eI z , with the identity operator 1 and < e <C 1. Throughout 
the paper we use the notation Ij = Oj/2 where a a are the familiar Pauli spin matrices. 

Radio frequency pulses are applied at the Larmor frequency 400 MHz of the protons. Experiments were performed 
with a Varian Infinity Plus NMR spectrometer applying a cycle time of t c = 2r = 100/iS and with relaxation time 
T2 ~ 200ms. The pulse sequence shown in Fig. [3] is based on the Carr, Pureed, Meiboom and Gill (CPMG)- 
sequence|22, 23] which was proposed almost 50 years ago for measuring Ti relaxation times in inhomogeneous fields. 
For our purpose we have modified this sequence such that the resulting evolution of the proton spins expressed by a 
sequence of echos which leads to the desired Gauss sum. 

The pulse sequence consists of an initiating 7r/2-pulse in y-direction which creates the initial density matrix p m = 
1/2 — e I x . This pulse is followed by a series of M + 1 7r-pulses with separation 2r which are individually phase shifted 
with respect to the a;-axis of the rotating frame by an angle |24l |25| . In the laboratory system the corresponding 
Hamiltonian reads 

M 

H(t) = hujQ I z + 2Trh^S(t - Tfe) cos(wr fc - (j>k) Ix (3) 

fe=0 

where loq and u denote the frequencies of the transition and the driving field, respectively. The pulses act at times 
Tfc = (2k + 1)t and the phases <pk will be chosen later in order to obtain a Gauss sum. 
With the help of the identity 

Uj(a) = cxp(-ialj) = cos(a/2) 1 — 2i sin(a/2) Ij (4) 

the Hamiltonian H within the rotating wave approximation, that is in the frame rotating with the frequency Ato = 
luq — lu takes the form 

M 

H = hAu I z + nh d(t — Tfe) (cos 4>k I x + sin <f>k I y ) (5) 

fc=0 

and the resulting time evolution operator 

U k = U,(AuT)U g {4> k )U x (ir)Ul(4 k )U,(AuT) (6) 

of the fc-th cycle reflects the three stages shown in Fig- 01 (*) Free time evolution for the time r followed by (ii) a 
7r-pulse which imprints the phase and (Hi) another free time evolution during the time r. The 7r-pulse eliminates 
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the effects of the free time evolution and ensures that the decoherence due to an inhomogeneous distribution of local 
fields is compensated by refocussing. 

At times r m + t we measure the polarization 

Tr (I xPm ) 
Tr(4p ln ) 

in x-direction with p m = p(r m + r) = UmpiJA^ where U m = U (r m + r) = JlfeLo ^fc ^ s the product of the time evolution 
operators Uk due to all previous pulses. 
When we apply Eq. Q to Eq. JSJ we find 

^ = H)Lil n (8) 
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which yields for the sum ,§( M ) = X)m=o Sm over ^ ne signals s m given by Eq. J7J) the expression 

M / m \ 
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Hence, the spin dynamics expressed by the signal s m depends on its complete history, that is the phases of all previous 
pulses. In particular, they enter as an alternating sum. It is here that the sum Eq. of the young Gauss comes into 
play. When we compare S^ M ' to the Gauss sum of Eq. J2J) we recognize that for the choice 
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for k = 



of the phases 4>k Eq. © takes the form 



with4 M) =Re^ M) . 

In Fig. 0] we display the results of our NMR implementation of our factorization scheme based on Gauss sums 
for N — 157573 = 13 ■ 17 ■ 23 • 31. On the top we show the time evolution of the spin under the influence of the 
particular pulse sequence given by the Hamiltonian H and the phases <j>k defined by Eqs. © and (fTtHl . respectively. 
As a measure of the dynamics we show the echo signal s m following from Eq. JJJ). For factors of N such as £ = 17 the 
signal is constant. Consequently we find for the average C^\i = 17) a value close to one as indicated in the bottom. 
However, for non-factors such as £ — 18 the echo signal oscillates around and leads to a rather small average value 
C [ ™\l = 18), shown by the arrow. We emphasize that due to the quasi-random interference of Gauss sums M = 11 
terms arc sufficient to discriminate factors from non- factors. 

Although decoherence is not a limiting factor in our scheme it is still interesting to investigate its influence. We take 
incoherent processes into account phenomenologically by introducing a Ti relaxation process and Eq. 10 is modified 
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w = mtt £ exp r m l) cos ( w t) • (12 
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We find that even for an appreciable decay, that is for 2Mt/T2 = 2 at the end of the sequence the pattern shown in 
Fig-fflby circles looks almost identical to the one without decoherence indicated by black dots, except for the reduced 
scale of the signal strength. This result is surprising since the signal sm at the end of the sequence has decayed to 
13.5% of its initial value. Hence, decoherence does not significantly influence our ability to distinguish factors from 
non- factors. 

In summary we have presented the first experimental implementation of a factorizing algorithm based on Gauss sums 
with an ensemble of single qubits. We have exemplified this technique factoring the number N = 157573. However, 
we claim that extensions to larger numbers are readily possible since our method capitalizes on the quasi-randomness 
of the phases of Gauss sums. As an outlook we factor in Fig. [3] the 24-digit number 



N = 1062885837863046188098307 
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FIG. 4: Experimental realization of factoring N = 157573 = 13 ■ 17 ■ 23 • 31 using the NMR implementation of the Gauss sum 
the top and bottom we depict the echo height s m measured 
respectively. For factors such as I = 17 the signals a 
average value Cjy (17) close to unity. In contrast, for a non-factor such as 
over m almost averages out as indicated by the arrow 
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measured at times r m + r and the resulting average C^ 1 - 1 for different 
are approximately constant as a function of m with an 
18 s m oscillates around zero and when summed 



with only M — 200 pulses. Since in NMR it is possible to even have up to 10 4 pulses within the decay time Ti several 
new questions emerge: (i) What is the optimal number M of terms in the Gauss sum, that is how many pulses are 
needed for a given N in order to discriminate factors from non- factors? (ii) How to overcome pulse errors? and (in) 
how to employ entanglement in order to achieve a speed-up? 

Obviously, the answers require a more detailed analysis. Therefore we can only speculate. However, our numerical 
experiments suggest a logarithmic dependence of M on N. In addition, pulse error correction techniques based on 
optimal control theory offer themselves. Finally, the entanglement of two and more spins might open a possibility to 
reduce the complexity. 
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FIG. 5: Probing the limits of our factorization scheme. Factorization interference pattern for the 24-digit number N = 
1062885837863046188098307 = pq obtained from a numerical simulation of the Gauss sum in the immediate neighborhood 

of the prime factor p = 790645490053 for M = 200. 
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by the Gauss reciprocity law0, Q. 



